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f^ , In this paper we present a new derivation of the QCD factorization. We deduce the fcr- and 

^Nj ' colhnear factorizations for the DIS structure functions by consecutive reductions of a more general 

theoretical construction. We begin by studying the amplitude of the forward Compton scattering 

IjT^, off a hadron target, representing this amplitude as a set of convolutions of two blobs connected 

ry\ • by the simplest, two-parton intermediate states. Each blob in the convolutions can contain both 

the perturbative and non-perturbative contributions. We formulate conditions for separating the 

CO , perturbative and non-perturbative contributions and attributing them to the different blobs. After 

that the convolutions correspond to the QCD factorization. Then we reduce this totally unintegrated 

(basic) factorization first to the kx- factorization and finally to the coUinear factorization. In order 

' (-H ' \ to yield a finite expression for the Compton amplitude, the integration over the loop momentum 

f^ I in the basic factorization must be free of both ultraviolet and infrared singularities. This obvious 

I ' mathematical requirement leads to theoretical restrictions on the non-perturbative contributions 

Oh' (parton distributions) to the Compton amplitude and the DIS structure functions related to the 

Qj , Compton amplitude through the Optical theorem. In particular, our analysis excludes the use of 

the singular factors a;"" (with a > 0) in the fits for the quark and gluon distributions because such 

factors contradict to the integrability of the basic convolutions for the Compton amplitude. This 

rf\ . restriction is valid for all DIS structure functions in the framework of both the kr- factorization 

^ ■ and the colhnear factorization if we attribute the perturbative contributions only to the upper blob. 

^\ ' The restrictions on the non-perturbative contributions obtained in the present paper can easily be 

^SJ , extended to other QCD processes where the factorization is exploited. 
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I. INTRODUCTION 

Factorization is the basic concept for providing the theoretical grounds for the use of perturbative QCD in the 
analysis of hadronic reactions. According to it, QCD -calculations can be done in two steps: first, the perturbative 
effects are accounted for with the calculations of Feynman graphs involved and second, the perturbative results are 
complemented by the appropriate non-perturbative contributions. There are two basic kinds of the factorization 
available in the literature: the collinear factorizationjlj and kx -factorization introduced in Ref. [3]. In the framework 
of the collinear factorization, the border between those stages is established by introducing a factorization scale /j, 
which at the same time acts as the infrared (IR) cut-off to regulate the IR divergences of the perturbative contributions. 
This kind of factorization is used in DGLAP[3| where the transverse parton momenta ki ± in the perturbative region 
are ordered Q as: 

n'^ < kl j^ < kl ^ < ... < Q^ . (1) 

Due to the DGLAP -ordering ([Ij, the virtualities of the initial partons are smaller than the virtualities of the 
ladder partons. The ordering ([iJ and the DGLAP evolution equations were introduced in the kinematics where 
X = Q^/2pq < 1. 

In order to obtain a generalization of the DGLAP equations in the small- a; region, the ordering ([T|) for the virtual 
quark and gluons in the perturbative region should be replaced by the ordering obtained in Ref. [^] : 

I? < -^ < -^ < ... <w = 2pq . (2) 

Pi P2 

We have used in Eqs. ([Tl [2]) the standard notations: p stands for the moment of the initial hadron, q is the moment 
of the virtual photon, Q^ = —q^, Pr {r = 1,2,..) stand for the fractional longitudinal momenta of the ladder partons 
and the factorization scale fi^ is the starting point of the Q^ -evolution. The ordering ([1]) is used for accounting for 
the leading logarithms of Q^ to all orders in as while ^ is used for the total resummation of all leading logarithms 
regardless of their arguments. The difference between the orderings H]) and ^ leads also to different treatments of 
as (see Refs. JL5j-[19] for detail). 

The common feature of Eqs. (H]) and (HJ is that the transverse momenta ki ± are restricted from below. This is 
necessary in order to regulate the infrared (IR) divergences arising from the double-logarithmic (DL) contributions. 
In this case /i acts also as an IR cut-off. Evolving the scattering amplitudes with respect to /i is the essence of the 
Infrared Evolution Equations (IREE) obtained first in Ref. [5[ for the quark scattering and then generalized to various 
other high-energy processes (see e.g. Ref. Q). In particular, a generalization of DGLAP was obtained in order to 
describe the DIS structure functions gi and F^^ at arbitrary x and Q^ (see the overview of these results in Ref. [3] ) . 

DL contributions are absent in the BFKL equations^, so this approach is IR-stable. As a consequence, the 
transverse momenta of the virtual partons here can be arbitrary small. On the other hand, the initial and final 
partons (gluons) in this approach are essentially off-shell. So, instead of /i, the transverse momenta k± of the external 
gluons act in BFKL as a new factorization scale, i.e. the border separating the perturbative part from the non- 
perturbative one. Such a factorization (the kr -factorization) was suggested in Ref. [2| . The value of k± is arbitrary, 
so the kT -factorization involves the integration over k± . 

The value of the IR cut-off fi is arbitrary (except for the requirement ^qcd < A*^ ^ Q^) but the final expressions 
for the scattering amplitudes and structure functions should be insensitive to /U, though both perturbative and non- 
perturbative contributions taken alone depend on it. The sensitivity of various physical observables to IR cut-offs has 
been a subject of a great interest among theorists. However, such investigations were often focused on the perturbative 
parts of the observables, leaving the non-perturbative parts unconsidered (see e.g. the recent overview Q and Refs. 
therein). In contrast, we consider in the present paper the IR- properties of both perturbative and non-perturbative 
parts. 

As it is well-known, the factorization means that any scattering amplitude and any DIS structure function can be 
represented as a convolution of the perturbative and non-perturbative blobs. In the present paper we consider the 
photon- hadron (Compton) scattering and begin with analysis of the forward (at i = 0) Compton amplitude A^,^. 
Representing A^^, as a convolution of two blobs depicted in Fig. [1] we study its integrability and then, using the 
Optical theorem, we apply the obtained results to the DIS structure functions. In contrast, most of the preceding 
approaches (see e.g. Ref. [13]) addressed directly to the DIS hadronic tensor W^i/ and because of that they could not 
obtain our results and arrive at the conclusions we make in the present paper. The intermediate particles in Fig. [T] 



^ The numeration of fcj x in JT] [2]l runs from the bottom to the top of the perturbative ladders. 



can be quarks or gluons. The number of them can be arbitrary. We consider only two-particle states connecting the 
blobs. On one hand, it is the simplest option; on the other hand, it corresponds (at large Q^) to the leading twist 
approximation. All blobs in Fig.[l]are not cut. The upper blob in Fig.[l]is supposedly obtained with any perturbative 
approach, including DGLAP or BFKL, or IREE evolution equations, though strictly speaking it can contain also 
unperturbative contributions. We study this issue in detail in Sects. II, III. 
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FIG. 1. Representation of A^i, through the convolution of the perturbative and non-perturbative blobs. 

The lower blob T is also a mixture of the perturbative and non-perturbative contributions. When the radiative 
corrections to the upper blob are neglected and the blob is considered in the Born approximation (see Fig. [5]), the 
lower blob can be regarded as altogether non-perturbative. For example, this takes place in the parton model. With 
the QCD radiative corrections taken into account, both the upper and lower blobs in Fig. [T] acquire the perturbative 
contributions. Nevertheless, following the standard terminology, throughout the present paper we address the lower 
(upper) blob as the non-perturbative (perturbative) blob. When the graph in Fig. [1] is cut in the s -channel, i.e. 
when all s- channel intermediate particles are on-shell, the lower blob = SsT corresponds to the probability to 
find a constituent (a quark or a gluon) in the incoming hadron. Throughout the paper we will address T and ST 
as the unintegrated parton distributions. When the graph is not cut, we define T as the result of applying the 
Dispersion Relations to ST. In the present paper we use the Feynman gauge for virtual gluons. On the other hand, 
we would like to stress that we consider the inclusive processes only, which makes us free of the analysis of the rapidity 
divergencies arising when the semi-inclusive processes are investigated. This kind of divergences was first investigated 
in Refs. [Ill [13 ■ See also recent papers [1^ and references therein. 

The available in the literature expressions for the unintegrated parton distributions were obtained without using 
any theoretical grounds. The only argumentation for fixing them was to fit the experimental data. In the present 
paper we obtain certain theoretical restrictions on the distributions following from the integrability of the convolution 
in Fig. [TJ i.e. from the obvious mathematical requirement that the integration over momentum k must be convergent. 
In particular, we show that such restrictions are compatible with the use of the singular factors ^ a;^° in the fits 
for the initial parton densities for the singlet structure function i^i, providing a < 1, and exclude such factors in the 
initial parton densities for the other structure functions. 

Our paper is organized as follows: In Sect. II we consider the convolution shown in Fig. [T] It represents the forward 
Compton amplitude A^i, in a factorized form. We express this amplitude in terms of the invariant amplitudes and, 
using the Optical theorem, we relate them in the standard way to the DIS structure functions. In order to make our 
reasoning simpler, we discuss in Sect. Ill the integrability of the convolution in Fig. [T] when the Born approximation 
is used for the upper blob. In Sect. IV we consider the convolution in Fig. [1] beyond the Born approximation. In this 
Sect, we reduce this convolution to the convolution corresponding to the QCD factorization where the upper blob 
is free of non-perturbative contributions. This factorization involves convolutions with totally unintegrated parton 
distributions, so we name it the basic form of factorization and discuss the similarity and difference between it and 
the standard (fcr and coUinear) factorizations. In Sect. V we study the impact of the radiative corrections on the 
integrability requirements obtained in Sect. III. In Sect. VI we show that the basic factorization can be reduced to 
the fcT-factorization only approximately and formulate the condition for such a reduction. In Sect. VII we continue 
studying the kx- factorization and consider restrictions on the parton distributions. In Sect. VIII we reduce the fc^- 
factorization to the coUinear factorization and formulate restrictions on the singular factors ^ a;^° in the DGLAP fits 
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for the initial parton densities for the DIS structure functions foUowing from the integrabihty of the basic convohitions. 
FinaUy, Sect. IX is for concluding remarks. 

II. REPRESENTATION OF THE FORWARD COMPTON AMPLITUDE AS A CONVOLUTION 

We start by considering the forward Compton scattering off a hadron. The amplitude Afj^i, for this process includes 
both the perturbative and non-perturbative contributions. According to factorization, Af^^, can be represented as the 
sum of the convolution shown in Fig. [T] where the intermediate partons are either quarks or gluons: 

(3) 
with 

MS{q,k)kT('lHk,p,m^„Sh) (4) 

Al%,iq,k)Tl,^Jik,p,m^,S^) 

where we have used the superscript q (g) to mark that the intermediate partons in Fig. [T]are quarks (gluons). In 

Eq. dH) Al^} and A^^jp^ stand for the upper blobs in Fig. [T] and r(«), T^P denote the lower blobs. Following the 
conventional notations, we have used the subscripts fi and v for the polarizations of the external virtual photon; p, a 
mark polarizations of the intermediate gluons; q,p and k are momenta of the photon, hadron and intermediate virtual 
partons (quarks or gluons) respectively, mh and Sh are the hadron mass and spin. In what follows we will focus on 
the integrations over fc, so we skip nifi and Sh in all subsequent formulae. Throughout the paper we will keep the 
standard notations Q^ — —q^, w — 2pq and x = Q^/w. 

Obviously, even at large Q^ and w the amplitudes Ajj',;, Ajfjpa can acquire unperturbative contributions at small 
k^. In order to keep these amplitudes completely perturbative, such a soft region should be excluded from Eq. ^. For 
example, in the framework of the collinear factorization it is done by introducing the factorization scale. In contrast, 
we will do it with imposing restrictions on the k^ -dependence of T^'^>,Tp^ but even before doing so we will address 
Ap'^ , Apjpa as perturbative objects. Besides them, the convolution in Eq. (U) involves the amplitudes T^''' and Tpa 
corresponding to the lower blob in Fig.[T]where non-perturbative contributions are collected. Strictly speaking, T^''^' 
can include perturbative contributions as well. However, these perturbative terms are quite similar to the perturbative 
contents of the amplitudes A/fj , Apjpcr . Throughout the paper we focus on the non-perturbative contents of T^^^ , Tpa , 
so we address them as a non-perturbative objects. As it is known, the perturbative amplitudes ApJ{q, k), Apjpa{q, k) 
can contain IR-divergent contributions. The virtual parton momentum k regulates these divergences and k^ (or 
k\) acts as an IR cut-off. We stress that the convolutions in Eq. @, though correspond to a factorization of Ap^ 
into the blobs, do not coincide with either the collinear or /ct- factorizations. Moreover, they do not imply any 
separation of the perturbative contributions from the non-perturbative ones, which is the fundamental concept of 
the QCD factorization. Because of that we will not apply the term "factorization" to the convolutions in Eq. ^ 
before obtaining conditions to separate the perturbative and non-perturbative contributions. Instead, we call them 
the primary convolutions. 

The imaginary part (with respect to s = {q + p)^) of Apy is related to the hadronic tensor Wp^: 

Wp, = -5sA^. . (5) 

IT 

In order to simplify the tensor structure of Ap^, one can use the conventional DIS projection operators PpJ enlisted 
in Appendix B. Using them, we represent Ap^ in terms of invariant amplitudes A^.: 

Apu^Y.Pl^'^^r. (6) 

Therefore every DIS structure function /^ is related to the invariant Compton amplitude A^: 

U = -^Ar. (7) 

TT 

Obviously, each invariant amplitude A^. can be expressed through the primary convolutions of the perturbative 
invariant amplitudes A)?'^ and the non-perturbative invariant amplitudes T*'^. Borrowing the terminology of the 



Regge theory, we can state that there are the singlet (As) and non-singlet (Ajvs) invariant amplitudes contributing 
to A^^. They have the vacuum and non- vacuum quantum numbers in the t- channel respectively. It is important 
that they exhibit different behavior with respect to s: As have an extra power of s compared to A^s- Generically, 
we can write the primary convolutions @ for them as follows: 



As = As®Ts, Ans^Ans^Tns- 



(8) 



Applying Eq. (jT]) to the singlet amplitude As yields the singlet structure function Fi whereas the non-singlet am- 
plitudes Ans lead to other structure functions. We would like to stress that this terminology is not quite adequate. 
Indeed, among these "non-singlets" there are the amplitudes with imaginary parts corresponding to the flavor singlet 
components of the structure functions (for example, the singlets (71,2)- 

III. COMPTON AMPLITUDE A^^ IN THE BORN APPROXIMATIONS 

In this Sect, we consider the primary convolutions ^ in the simplest case when the perturbative amplitudes 
Ap'u, A\fjpa in Eq. (g]) are calculated in the lowest-order approximation. The convolution for A^^, is depicted in 
Figs. 2,3. It is convenient to consider this approximation for A}^^ and Ajfjpa severally. 

A. Lowest-order approximation for Af^J 

The case when the amplitude A}^,! is calculated in the Born approximation is depicted in Fig. [5] 
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FIG. 2. Born approximation for the amplitude of the forward Compton scattering. 



It corresponds to the use of the parton model. Dropping the QED quark-photon coupling e^ as an unessential 



factor, we can write the Born perturbative contribution Aj^i, as follows 
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where we have neglected the quark mass. Therefore, kK^ly k can be written as follows: 

^" ik + q)^+ie [k-qf+ie 



(9) 



(10) 



The factors B^^^^q, k) in Eq. (llOp are considered in detail in Appendix A. Before performing a detailed calculation, 
let us give a simple Euclidean-like estimate for the ultraviolet behavior of T'^^ in Eq. ^ when Ajfj = is given by its 
Born value. Obviously, we can write d^k — k^^dkdfl^. Then, at large k the Born amplitude in Eq. © Ap'i, ^ 1/k^ 
and therefore at large k 
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This integral is convergent at large k only if 
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with h> . 

Now let us investigate this case more carefully. It is convenient to use the Sudakov variables introduced in Ref. |jU 
We will use them in the following formn: 



k — ~aq' + Pp + k±, 



with q' — q + xp, x = Q'^ /w, so that q'' 



and 



k = —wa(i — fc_L, 2pk = —wa, 2qk — w{l3 + xa), {q + k) — ■w{(3 — x){l — a) — fcj 

(r) 



(13) 
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B^i, can be simplified with using the projection operators P/lJ defined in Appendix B. Then we arrive at the 
following expressions for A^^: 



Ai'} = E ^il^^^? = E ^M^^ / dkldl3daAi'^ ^) (q, k) 
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where Ar are the invariant amplitudes in the Born approximation and B is defined in Eq. (jASp . We have also 
provided the non-perturbative amplitudes T' with the superscript r as they can be different for different ArT'^. 
Obviously, Tr {k,p) depend on the invariant energy and the external virtualities: 



T,(«)(fc,p) = T,(«)((fc +p)2, fc2) ^ T^q) ^^a, {wafi + kl)) 



(16) 



The Born invariant amplitudes A)? are well-known. For example, the amplitudes Af and ^3 related to 
the structure functions Fi and gi respectively are 
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Taking the imaginary part of Jv-^ , A3 with respect to w and neglecting the virtuality /c^ 



j'(peri i3) 



arrive at the well-known result of the parton model: S^i 3 ^ <5(/3 — a;). 

Now let us integrate Eq. (|15p with respect to a and focus on the region of large \o\. In this region the amplitudes 
A^ ^ ^ a/a, both the denominator and B in Eq. p^ are ~ a^, so the a-integration looks at large \a\ as follows: 
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This integral is convergent at large |a| only if 
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with h > 0. Eq. (TTOl) confirms the estimate made in Eq. p2)) . According to Eq. (TTS)) . the invariant energy s' = {p— fc)^ 

of T^*"^ at large a is s' = wa — waP — fc^ w wa. So the meaning of Eq. ([T9| is that Tr (s') should decrease faster than 

1/s' with growth of s' in order to prevent UV divergency in Eq. (jlSp . Obviously, AJ^' has non-vacuum quantum 
numbers in the t- channel and therefore it contributes to the non-singlet part of the Compton amplitude Af^^ . 



We account for Q but neglect the quark mass. 




» 



q q 



p p 



Q q 



q 




p p 



(b) 




c) 



FIG. 3. Gluon contribution to the forward Compton scattering in the lowest order. 



B. Lowest-order approximation for A 



(9) 



The case when the amplitude yljfi; is calculated in the lowest order is depicted in Fig. [31 

For the shortness reason, we will address it as the Born approximation as well. The primary convolution for the 
Compton amplitude A]^^ in this case is 
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(20) 



where we have used the notation T^^, pa-{q, k) for the upper, perturbative blob given by the quark boxes, one of them 
depicted in Fig. [31 After applying the projection operators, Eq. (|20p can generically be rewritten in terms of the 
invariant amplitudes as follows: 



A(3)(g,p) 



^^pa(q,k)-^ 



i2ny 



[k^ 



T^'Ak,p). 



^2 P<r 



(21) 



-(5) 



The rough estimate for Tpa can be obtained similarly to the one in Eq. (|12p . Writing d k — k dkdfl^ and assuming 
that k is large leads to 
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The integration over k is free of the UV singularity when 
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which differs from Eq. (|T2|) . 

Now let us make a more detailed estimate of the same Born case, using the Sudakov parametrization (|13p for k. 
At large a Eq. ([21]) is 



A(9) 



j daT{q,k)p,^T^/J{a,e). 



(24) 



The contribution a^ in Eq. (|24p comes from k^ . According to Eq. ([TS]) . fc^ = —wafi — kfj_, i.e. fc^ ^ a at large 
a. The perturbative factor Tp„ is known to be free of the UV singularities and can be calculated by integrating over 
momentum fci (see Fig. [3]). This integration yields different results for the singlet and non-singlet. Also a dependence 

of Tpa on a can be different for the singlet and non-singlet. 



In the singlet case the polarizations p and a are longitudinal and Tpa can be represented as follows: 



Tl^J ^'^-^T'i\a,k% (25) 

so the integration over a in Eq. (|2ip is reduced to 



This integral is convergent when 



^s - I da—T'g'>ia,k'). (26) 



T^f'> ^ a-'\ (27) 



It coincides with the rough integrability requirement in Eq. (|23p but differs from Eq. (J12p for the quark contribution. 

The non-singlet amplitudes Ajvs are either flavor non-singlets or spin-dependent. In the latter case Tpa is anti- 
symmetrical in pa and can be represented as 

T(,sJ = iIIllLep^xrSxkrT^^'lia, fc^), (28) 

with Sx, rrih being the spin of the hadron target and the mass scale respectively, or in a similar form. Therefore, in 
the non-singlet case Eq. ((27|) is replaced by 

da'^T^j^l(a,e). (29) 



»2 
a 

This integral is convergent when 



r^'s - «"'"'. (30) 

which coincides with the integrability requirement for the quark non-singlet case. Therefore, the integrability require- 
ments (|19l30p for the non-singlet part of A^^, in the lowest-order approximation do not depend on the kind of the 
intermediate partons. 

IV. COMPTON AMPLITUDES As, Ans BEYOND THE BORN APPROXIMATION 

When the perturbative amplitudes As, ^jvs in Eq. ^ are calculated in the Born approximation, the amplitudes 
T5, T/vs can be regarded as completely non-perturbative objects. Studying As, Af^s beyond the Born approximation 
brings the radiative corrections to As, Ajsis- In principle, the radiative corrections can also be placed into Ts, T]\js- 
It converts the Born invariant amplitudes A[^^ of Eq. (|17p and Tp„ of Eq. pop into more involved invariant am- 
plitudes which we generically denote ^5 ^5- This can have an impact on the a-dependence of the non-perturbative 
blobs Tgf^g and violate the integrability requirements for the singlet ([27l) and non-singlet ()19|30p amplitudes. In 
the previous Sect, we showed that the integrability requirements do not depend on the kind of the intermediate 
partons in the convolutions, so in the present Sect, we focus on the quark amplitudes A^ and will skip the su- 
perscripts q, g in the sequential formulae. Each of As^ns can depend on a through fc^ only. Besides, they include 
the infrared-sensitive radiative corrections becoming singular at small k^. In particular, there are the logarithmic 
terms ^ ln"(w^/A:^),ln((5^/fc^). In addition to them, the first-loop radiative corrections to the singlet unpolarized 
amplitudes As yield the infrared-sensitive power term w/S/k'^. As it is known, such a term is originated by the 2-gluon 
intermediate state, with the gluons having the longitudinal polarizations. Such amplitudes can be written as follows: 

As = (w/3/P)M5(ln(w/3/fc2), HQ^k'')). (31) 

As written explicitly, Ms in Eq. (1311) includes the infrared-sensitive logarithms. Of course, it can also include other, 
infrared- insensitive corrections. The non-singlet amplitudes ^^5 do not have the power term (w/3/fc^), so we write 
them in the following way: 

Ans = MNs{Hwf3/k^), HQ^k^)). (32) 



In these terms we can rewrite the primary convolutions in Eq. ^ as foUows: 

As ^ J *i |.0 (f ) M. (H.m% HQye» j^j;;^ 
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(34) 



We remind that each of the amphtudes Ms.ns and Ts.ns in the convolution p3l34p can contain both perturbative 
and non-perturbative contributions, so these convolutions (we have addressed them as the primary convolutions) 
cannot be identified with the factorization convolutions. 



A. Converting primary convolutions into factorization convolutions 

Ms, Mms in Eqs. p3l34[) can be entirely perturbative objects only if the region of small k^ is excluded from these 
convolutions. Besides, the integrals in (|33l34p should be IR- stable. These goals can be achieved when 



Tns - ik') 



2N|77 
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and 



Ts 
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2n1+»7 



(36) 



with ?7 > 0, at small k^ . 

Imposing conditions p5l35p makes the amplitudes Ms^ns (i-e. the upper blob in Fig. [1]) free of non-perturbative 
contributions and it also solves the problem of the mass singularities [2 1|. However the lower blobs, Ts^ns can include 
both the perturbative and non-perturbative contributions. When the perturbative blobs were considered in the Born 
approximation, the lower blobs were regarded as an absolutely non-perturbative objects. In this paper we advocate 
the point of view that beyond the Born approximation the lower blobs contain non-perturbative contributions only. 
Indeed, let us add a radiative correction (for instance, a ladder gluon propagator as shown in Fig. |4^) to the Born 
graph. 





p p p' ^ p 

(a) (b) 



FIG. 4. Radiative corrections to the Born amplitude. 



Obviously, it can be included into the upper blob. This procedure can be repeated each time when a new radiative 
correction is added to Figs. 2,3, except for the case when the new propagators connect the upper and lower blobs as 
shown in Fig. HJd. Such graphs involve the i-channel intermediate states with three or more partons. Analysis of them 
requires a special attention (see e.g. Ref. [23]). We do not consider them in the present paper. 
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Eqs. P3I34[) complemented by the IR-restrictions (l35l |36| rigorously correspond to the QCD factorization concept 
where the non-perturbative contributions should be attributed to Ts,Tjvs while the amplitudes Mg, M^g can be 
regarded as entirely perturbative objects. Nevertheless, it is clear that the form of the convolutions in these equations 
does not correspond to the form of either the kr- or the collinear factorization. Indeed, these factorizations involve 
the same objects (the perturbative and non-perturbative blobs) as in Eqs. (J33l34p . however with the smaller number 
of the integrations. Let us remind that the expressions with the fcT-factorization include two integrations only: over 
/3 and k±, while the collinear factorization involves the integration over the longitudinal variable /3 only. In contrast, 
Eqs. p3l34p additionally include the integration over a. In order to distinguish the convolutions in Eqs. P3I34I) from 
the convolutions with the kr- and collinear factorizations, we will call the totally unintegrated form of the factorization 
convolutions in Eqs. p3l34p the basic form of factorization. Below we show how to proceed from the basic form of 
factorization to the k^- and collinear factorizations but before it, let us write down the factorization expressions for 
the DIS structure functions. 

B. Basic form of factorization for the structure functions 

Applying Eq. ([7|) to the invariant Compton amplitudes p4p complemented by the IR-restrictions ([55l [M| , we arrive 
at the following expressions for the non-singlet structure functions in the basic form of factorization: 

Ins = j dkl'^daf^^l^'^ {\n{wPie)M{Q^ 1^)) ^—^-^^^s{wa,e) (37) 

where fj^g^ = {l/Tr)^AlNS stands for the perturbative contributions to the non-singlets and ilNs{'wa,k'^) = 
'^Tj^s{wa, fc^) are the totally unintegrated non-singlet parton distributions. Similarly, we obtain the following repre- 
sentation for the singlet structure function from ([55)) : 

fs = f dkl^daft'''^ (ln(u;/3/fc2),ln(QVfc^)) . ^^ , , . , ».(^a, fc^), (38) 

where the perturbative singlet contribution is 

^(peW) ^ ^^ ^i/^)^Ms. (39) 

and Qs{wa, fc^) = '^Ts{'wa, fc^) stands for the totally unintegrated singlet parton distributions. 

V. CONDITIONS OF INTEGRABILITY FOR NON-SINGLET COMPTON AMPLITUDES 

Let us consider the integrations in Eq. p4p for the non-singlet Compton amplitude. Generally, the order of the 
integrations in Eqs. (I33l34p can be arbitrary. It is convenient to integrate over a first. This is usually performed by 
using the Cauchy theorem. So, the integration line — oo < a < oo is complemented by a semi-circle in the upper or 
lower half of the complex a -plane and after that singularities of the integrand with in the a -plane (poles and cuts) 
should be found. The integral over the semi-circle Cr vanishes when the radius of the semi-circle tends to infinity, 
so the result is determined by either residues at the poles or integrals along the cuts of the integrand in the a-plane. 
Now let us consider the integral over a in Eq. (|34)) . For integration along Cr, where |a| -^ oo, the a^ -terms both in 
B{k) and in the denominator of Eq. (pS)) cannot be neglected, so the integral along Cr behaves as 

Q,3 

da^TNs{a,kl). (40) 

It vanishes at |a| — >■ oo if 

Tns - a^^^'' (41) 



at large a, with arbitrary h > . This restriction coincides with Eq. (j30p obtained in the Born approximation. 
Therefore, accounting for the radiative corrections to the non-singlet amplitudes does not change the a -dependence 
of T at large \a\ and has the same meaning: T^vs should decrease faster than 1/s' when the invariant energy s' = [p—kY 
grows. 
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A. Integrability of the singlet Compton amplitude 

The singlet invariant amplitude As is defined in Eq. ([33|. The main difference between Eqs. (l33l) and the non-singlet 
amplitudes of Eqs. (IM|) is the factor wf3/k'^ in Eq. ([55)) . Therefore, the integration over a in Eq. ([55]) leads to new 
requirements for the singlet amplitude T5. Using the same arguments as for T^vs, we conclude that the integration 
over a at lal -^ 00 looks now as follows: 



/■ 



a 



2 



This integral is convergent when 



da^Ts{a,k'). (42) 



Ts ^ a-\ (43) 



which coincides with the integrability requirement (|27p for the singlet amplitude in the lowest-order approximation. 
It follows from Eq. ()13|) . that wa = 2pk is at large a the invariant energy of the non-perturbative blobs Ts, ns, so 
Eq. (H51) . similarly to Eq. (|1T|) . predicts the dependence of Ts on the invariant energy s' = wa at large s': Ts decreases 
with s' faster than s'^'\ 

VI. REDUCING THE BASIC FACTORIZATION TO kr -FACTORIZATION 

In order to bring the convolutions in Eqs. (J33l34p to the form of the fc^ factorization, the integration over a should 
be performed. Obviously, this integration cannot be performed in the straightforward way without dealing with the 
perturbative parts of the convolutions because they depend on fc^ and fe^ depends on a through Eq. (|13l) . The only 
way to bring these convolutions to the form corresponding to the fcy -factorization is to impose the following restriction 
on the longitudinal component wa(3 of fc^ compared to its transverse component: 

wa/S < k]^ . (44) 

Eq. dm allows us to disentangle the a and /3 -dependence in Eqs. P3l34l37l38p . Under this condition, the upper 
(perturbative) blob in Fig. [1] approximately depends on /3 and k± whereas the lower blob approximately depends on 
a and k±. It allows us to perform the integration over a in Eqs. p3|34|37|38|) without involving the perturbative 
components, which makes possible to bring these equations to the form corresponding to the fcy -factorization. To our 
knowledge, Eq. (j44p is not written explicitly in the literature. However, we would like to stress that it has commonly 
been used, though in an inexplicit way and in other term^j, in order to obtain integral representations of the DIS 
structure functions, including DGLAP. At a; ^ 1, the strong inequality in Eq. pi)) should be replaced by the weaker 
condition wa/S < k^. Eq. (HH) perfectly agrees with the well-known fact (but does not follow from it) that the 
virtualities of the ladder partons are transverse. Applying Eq. p4|) to A^s allows us to neglect the a -dependence in 
the infrared-sensitive terms in M^vs and to keep it in Tpfs only. Performing the integration over a in Eq. (I34p . we 
arrive at 

ANs{q,p) = y -^j jMms (ln(«;/3/fcl),ln(gVfci)) -^Nsiwli.kl), (45) 

with 

, „, „ /■°° BTNs(wa,wal3 + k^A fi'l/^P .9. 

^NS {wP, fci = fci / da ^, ^ ' ,,,, ^' « / daTNsiwa, kl). 46 

At a; '^ 1, the upper limit w of the k\ -integration in Eq. (|45l) can be replaced by Q^ as it takes place in DGLAP. 



For example, in the book | 22|) the Sudakov variables a and A;^ are replaced by 'm? = {p — k)^ and k^ respectively. 
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Taking the imaginary part converts Eq. (|45l) into the following approximate expression for the non-singlet structure 
functions /ns'- 

fNsix, Q') ^ r^ f' ^AT^ {w(3, Q^ fci) -s^^s {w(3, fci) , (47) 

"'0 «^_L J/3o P 



with 



and 



Pq^ X + k\/wKin\Six[x,k\/w\ (48) 



'i^Ns{wP,k\,m) p=i I da'^TNs{woi.,k\,m) = ^Ns{wl3,k\) (49) 

Jk\/w 

where we have introduced the non-singlet unintegrated parton distributions $Ars(/3, fc^). They correspond to the 
initial parton densities in the case of the k±_ -factorization. Now we can represent the non-singlet structure functions 
in terms of convolutions of the perturbative (partonic) and non-perturbative parts: 

fNs{x, Q') - T'^ f fAT^ {wP, g^ fci) $^5 {WP, fci) . (50) 

Jo '^± Jpo P 

The structure of Eqs. P51 [50]) is quite similar to the one in Ref. [3| and corresponds to the fcr -factorization. 
Introducing xq = k\/w and a mass scale ra (for example, it can be the hadron target mass) allows us to write 
Eq. ([50)1 in the following symmetrical form: 

fNs{x, Q') = r^ f' ^/(,^-*) {x/l3, gVfci) $^5 {xo/13, fci/m^) . (51) 

The attractive feature of Eq. ([5T|) is that it exhibits the remarkable symmetry between the arguments of fj^g^ and 
^NS- a; O Xo, Q^ -H> fci and fci ^ w?. Nevertheless, let us remind that though x and Q^ are conventionally regarded 
as independent variables, strictly speaking they are not independent. Indeed, they can be regarded as independent 
only in the case when Q^ is kept fixed while w is scanned. In the opposite case when w is fixed and Q^ is varied, 
the variables Q^ and x are not independent at all. So, the parametrization of the structure functions in terms of the 
really independent variables w and Q^ would be preferable at least for the theoretical analysis. 

Similarly to the non-singlet case, the singlet amplitude As can be written as follows: 

Asiq.p) ^£^Jj (^) Ms (ln(u;/3/fci), In(QVfci)) ^siwP, fci), (52) 

with 

9N . 9x9 Z""" BTs(wa,wal3 + k'^.) fkl/^f^ . o. 

M/S {wp,kl) = (fci 2 / da 7 ' /; ^' « / daTswa,kl). 53 

The unintegrated parton distribution ^{w/S, fci) can be defined similarly to the non-singlet case: 

fkl/wp 

9*s (w/3, fci, m) « / daSTs (wa, fci, m) = $s(w/3, fci). (54) 

Jk\/w 

Therefore, the singlet structure function fs (i.e. the DIS structure function Ff) is 

/s(x, Q') ^ r^ r fft"^'^ {wl3, g2, fci) $5 {yjl3, fci) . (55) 



13 
VII. INFRARED BEHAVIOR OF PARTON DISTRIBUTIONS IN THE fcr-FACTORIZATION 

The generalized parton distributions Ts^ns (in the expressions the Compton amphtudes) and ^s,ns (in the expres- 
sions the structure functions) include unperturbative contributions. When the kr -factorization is used, there is no 
problem about UV divergences but the IR divergences at small k± must be regulated. Below we study the impact of 
the integrability conditions on the properties of the parton distributions in the framework of the kx -factorization. 

Let us first consider the fcj_ behavior of the non-singlets. The explicit expressions for Mats' can be different, 
depending on the approach chosen to calculate it but their IR-divergent contributions in Afjvs are always logarithmic: 

Mjvs ^ In" («^/3/fci), In(QVfci), (56) 

with n,l = 1,2,... So, fc^ <C w/3, Q^ in any perturbative approach. It means that the k± -integration is free of 
ultraviolet divergences. However, the infrared divergences can appear in Eqs. (j45|50|) . In order to avoid them, we 
have to put a restriction on Tms at small k±: 

Tns - (klr (57) 

where 77 > 0. The main difference between Eqs. (|52I55I) and the non-singlets (I45l50p is the factor wf3/k'^ in Eqs. (|52l55p . 
It changes the non-singlet restriction (|57)) for 



Ts ^ (^1)^+" (58) 

at small k±. Besides regulating the IR-divergences, Eqs. (|57l58p have another important meaning for the kr factor- 
ization: they allow to neglect the region of small k± in the convolutions, solving therefore the problem of the mass 
singularities |2l|. 

The initial parton densities $ contribute to the DIS structure functions when the k± -factorization is used. In this 
case, k'j_/f3 is their invariant energy, so Eqs. (HTll43p . with a replaced by k\/wf3, predict the following dependence of 
the initial parton distributions on j3: 

^NS -^ /3\ (59) 

$s - p-^+'\ (60) 

They mean that ^ns should be regular in (3 when /3 — > and $5 can be singular but with h < \. These restrictions 
and the ones in Eqs. (|57l [58| must be respected when fits for the initial parton distributions $5, ^ns are composed 
in the framework of the fc^ -factorization. 

VIII. COLLINEAR FACTORIZATION FOR THE STRUCTURE FUNCTIONS 

In this Sect, we consider the infrared properties of the DIS structure functions when the collinear factorization 
is introduced. Historically, the collinear factorization 1] was the first one considered in the literature. In general, 
the basic assumption of the collinear factorization is to introduce the scale fi for the parton virtualities fcf , with 
kf > /i^ > Ag^^ as an explicit border between the perturbative (hard) and non-perturbative (soft) domains of 
QCD. For describing the DIS structure functions, this factorization first was used to construct the DGLAP evolution 
equations where the ordering ([1]) was suggested to evolve the structure functions from /i^ to Q^ . Later, the collinear 
factorization combined with the ordering ([2]) was used to describe some of the structure functions in the small-x region 
where the total resummation of In (1/a;) becomes essential (see Ref. |7| and Refs therein for detail). 

A. Transition from the kr -factorization to the collinear factorization 

Eqs. (|45I52I) for the forward Compton amplitudes as well as Eqs. (|50l55p for the structure functions involve integra- 
tions over j3 and fc_L, whereas similar expressions in the collinear factorization operate with integrations over (3 only. 
So, our aim now is to integrate over k± in Eqs. (|45I52I50I55I) without dealing with the perturbative contributions. 
However, both the perturbative and unperturbative factors in Eqs. (I45I52I50I55P depend on fc_L. It makes impossible 
any straightforward fcj^ -integration of the unperturbative terms. Therefore, a transition from the k±^ -factorization 
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to the coUinear one can be done only approximately. In a sense, it is similar to the transition form the basic form of 
the factorization to the k^ -factorization we have considered in Sect. VI. 

In order to do it, let us first suppose that ^s,ns{= ^) in Eqs. (|50I55|) can be represented as 



$ = ^(/3,fci)eo.2_fc^), 



(61) 



with ip{/3,k'^) decreasing rapidly at large k± and respecting the restrictions in Eqs. (|57l58p . This choice of $ is 
depicted in Fig. [SJ 



$(/3,fcl 



FIG. 5. The plot for $(/3,ifcl) defined according to Eq. (jeT 



It means that fi is the strict border between the contributions forming the upper and lower blobs in Fig. [T] the 
perturbative corrections with k± > fx are attributed to the upper blob and participate in the orderings of Eqs. ([1] [5]) 
while the lower blob includes the contributions with k± < fi. 

Alternatively, the transition from the fc-p- factorization to the collinear factorization can be done without using the 
0-function. Let us assume (see Fig. [6|) that $ has a sharp maximum at fc^ = /i^ ^ Q^ saturating the integration 
over fc_L in Eqs. (|50l5ip . 

$(/3,4) 




FIG. 6. The plot for <E>(/3,A:i) defined according to Eq. (I62j 



For example, 
peaked factor). 



it takes place when the fc^- dependence of $ includes the Gaussian exponential (or a similar sharply 
$(/3, fc^^) can also include a flatter function of k\: 



$ = (/?(/?, fc^)exp -o-(fc^-^ 



(62) 



where the function (/)(/3, k"^) in Eq. (p^ respects Eqs. (|57l58p . In this case fi cannot be regarded as the border between 
the contributions to the upper and lower blobs in Fig.[TJ $ in Eq. (1621) comprises the contributions with arbitrary k±. 
Substituting Eq. ((6T|) or (|62l) in Eqs. ()45|52|50|55p allows us to integrate $ over fc^ independently of j(p'=''*)^ arriving 
at the standard expressions with the collinear factorization. In particular, combining ((6T|) or (|62|) with Eqs. ([51] [55]) 
we arise at the following generic expression for the structure functions: 
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(63) 



{pert) 



where /3o(/i) — x + ^P /w ^ x. We have used m Eq. ([M]) the generic notations /, fr both for the singlets and 

non-singiets; the subscripts q, g in Eq. (j63p refer to the quark and gluon respectively. The perturbative contributions 
fr^ in Eq. (p5|) are usually calculated either in the framework of DGLAP or with total resummation of the leading 

The transition from the /cT-fa-ctorization to the collinear factorization does not affect 
so the singlet (/5s(/3,/U^) and non-singlet ipNs{/3,fJ''^) parton distributions used in 



(64) 



(and sub-leading) logarithms 
the /3-dependence in Eqs 



the framework of the collinear factorization should obey the following relations respectively: 



^s(/3,/^')'-/3- 



-l+h 



(65) 



with h > 0. We would like to stress here that the restrictions (IMl 1551) are not related to the integrability of the 
collinear factorization convolution in Eq. (|63p in itself. Indeed, the integral in Eq. (1551) is convergent at any (p because 
the integration over /3 here runs over the compact region [(3o{iJ,); 1]. However, violations of (|M1[55)) would destroy the 
integrability of the more general factorization convolutions in Eqs. P3l34p . We remind that in the present paper we 
have derived ([Ml US]) from Eqs. p3l34l) . 

Let us notice that /i^ in Eq. ([551) stands both for the factorization scale and for the starting point of the perturbative 
Q^-evolution. On the other hand, this evolution can start at an arbitrary point. When it starts at Q^ = /i§ (with 
/iQ > /i^), Eq. ([63]) can be written as 



fix,Q' 



r=q;g 



Mp. 



^fi^^^'^x/(3,Qy^^l)^^^\p,^^l) 



(66) 



with (p^'^'{J3, ^Iq) being new parton distributions. They are defined at the new factorization scale /ig. The values 
of Hq are located within the range of the Perturbative QCD, so the parton distributions (^*^''^(/3, /ip) include both 



{pert) 



non-perturbative and perturbative contributions. Although both f!f^'^^' (x//3, (3^//ig) and (^'^'''(/3, /Iq) depend on /ig, 
the integral in Eq. (|66p is obviously /ip-independent exactly as it is independent of contributions of any intermediate 
integration point. Indeed, the inclusion of the perturbative contributions into (f>{x, /ip) implies that there is another 
scale /i (with /i < /iq) in ^(a;, /ig), so the perturbative methods can be used for the evolution from /i^ to /ig. Obviously, 
the transition from (fis.Nsix, /i^) to ips^Nsi^j Mo) can be done with the same technique as the evolution from /ig to Q^ 
in fipt^rt). 



/(Q') = uiQ', fil)^ ® ii^l), ^{til) = u{^Jil m') ® ^ip^)- 



(67) 



where U{b, a) is a generic notation for the evolution operator, with a{b) standing for the starting (final) point of the 
evolution. We have skipped dependence on all variables in Eq. (p7)) inessential at the moment. For instance, in the 
case of the LO DGLAP ^(/^g) ^^'^ vil^'^) a-re related in the momentum space as follows: 



v{^,tJ'l) = 



ln(/i2/A2) 



-|7<'>M/b 



v{^,fJ''^), 



(68) 



where 7'-^-'(w) is the LO anomalous dimension and b = [33 — 2n/]/(127r) is the first coefficient of the /3 -function. In 
this case the perturbative structure functions f^P'^^*\uj,Q^) and /^^'^'"'-'(aj, /ig) are related similarly: 



expert) 



(aj,( 



"in (gVA2) 

ln(/.2/A2) 



7(i'(<^)/& 



/(P-*)(c.,/i2) 



(69) 
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Substituting Eqs. (J68l69p into (|66|) kills the dependence of f{x, Q"^) on the arbitrary scale /Iq but does not affect 
its dependence on the scale /i. We remind we have defined the scale jjl in Eqs. (|61l W^ as the scale separating the 
perturbative and non-perturbative contributions. 

Our restrictions (|64l65p look invalid in the framework of the approach suggested in Ref. [23]. We will address 
it as the EGR-approach. In this approach the factorization scale /lo can also be chosen arbitrary and the QCD 
coupling is kept fixed. The EGR-approach is based on the observation that the lader Feynman graphs contributing 
to the Gompton amplitude yield the sub-leading IR-stable logarithmic contributions in the kinematic region where 
the virtualities fc^ (r = 1, 2, ..) of the vertical partons are very small: 

Q<\kl\<^ll (70) 

while virtualities of the horizontal partons are greater than /i§. Like in Eqs. ()1|2|) . we have numerated the ladder 
momenta kr in (|70|) from the bottom to the top of the ladders graphs. Integrations of the ladder graphs over kr in the 
region (j70p yield logarithmic contributions ^ In" x for the non-singlets and ^ (I/2:) In" x for the singlet. Resummation 
of them leads to the Regge factors ^ a;^° for the non-singlets and ^ a;^^~° for the singlet (see Refs. [23,[2J| for detail). 
As fc^ in the region ([70]) are small, these Regge factors, according to Ref. [23] should be included into (^(x, /Iq) . The 
sub-leading logarithmic contributions ~ In" x are IR-stable, so the EGR-approach operates with the scale ^0 only 
and does not need any additional scale /z. This approach obviously contradicts our approach and makes invalid our 
restrictions (|64l65p . We plan to compare our and EGR approaches in detail in our forthcoming paper. Here we just 
mention that the EGR- approach leads to a much more complicated factorization construction than the one depicted 
in Fig. 1: instead of two blobs it involves a pile (etagere) of perturbative and non-perturbative blobs. 

B. Restrictions on the singular factors in the DGLAP fits 



Now let us focus on the convolution in Eq. (163)) . with the DGLAP expression used for the perturbative structure 
functions fr^ ■ In the standard notations, Eqs. (1551) takes the following form: 

/(:., Q')=Y. tl7 f DGLAP i^/y, 0"!^^) Sr{y, A.^) . (71) 



r=q,g 
,2 



where we have dropped the term 11 /w in the lowest integration limit. The initial parton densities 6q{x),6g{x) in 
Eq. (j71l) are obtained by fitting the experimental data at /i ^ 1 GeV. They include both singular factors a;"" (with 
a > 0) and regular factors ~ (1 — x)^, (1 -I- cx"^), with b,c,d> 0. 

Comparing Eqs. (|63l) and (|64l65p with ([7T|) drives us to conclude that the standard DGLAP fits with the singular 
factors x~°' for the non-singlet structure functions are excluded by Eq. (IMl) . On the other hand, Eq. (I55|) admits the 
use of the singular factors for the singlet Fi , however with the exponent a < 1. We would like to stress here that the 
use of the singular factors cannot be excluded by requirements on the integrability of the structure functions (where 
all integration limits are finite). However, the use of such factors contradicts to the integrability the expressions for 
the forward Gompton amplitudes where the integrations run over the whole phase space. Now let us comment on 
these results in more detail and consider below both the fcy- and collinear factorizations. 

First of all, we infer that the fits for the polarized quark and gluon distributions should satisfy the restriction ([59]) 
or (|64p , depending on whether the kx- or collinear factorization is used. 

Then, the spin-independent non-singlet structure functions involve the unpolarized quark distributions only, so each 
of the quark distributions should again satisfy Eqs. (|59l64p . It is true also for the singlet F2. 

The difference between the restrictions (|59l64p for the non-singlets and the ones (|60l65p for the singlets arises 
entirely from the perturbative components of the singlet Fi (= Ff) and the other structure functions. Ff involves 
the unpolarized quark and gluon distributions. The quark distributions in Ff are the same as the ones to the non- 
singlets Fi, so they satisfy Eqs. (I59l64p but the restriction on the gluon distribution to Ff is softer because it is 
given by Eqs. (|64l65p . On the other hand, the same gluon distribution is supposed to contribute to the singlet F2 and 
therefore it should also satisfy Eqs. ([5^[M|) . Eventually, we conclude that the singular factors x^" should be absent 
in all parton distributions. 

In Refs. [3, [20] we showed that the singular factors in the DGLAP -fits were introduced to match experimental data 
at small x. The reason for introducing them is that they cause a steep rise of the structure functions at small x, thereby 
making possible the use of DGLAP in the small- a; region. As a matter of fact, they mimic the total resummation of 
the leading logarithms of x, which is beyond the reach of DGLAP. With the resummation taken into account, such 
a rise of the structure functions at small x is achieved automatically, so the singular factors become unnecessary. In 
addition, we have shown in the present paper that there are theoretical restrictions on the singular factors following 
from the requirement of integrability of the factorization convolutions for the forward Gompton amplitudes. 
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IX. SUMMARY 

In the present paper we have derived the kx- and colhnear factorizations for the Compton amphtude and DIS 
structure functions and obtained the restrictions on the parton distributions to the DIS structure functions. The 
parton distributions available in the literature are fixed phenomenologically. In contrast, we have obtained the 
theoretical restrictions on them, exploiting the obvious mathematical requirement of integrability of the factorization 
convolutions. 

We began by considering the convolutions ^ depicted in Fig 1 where the amplitude of the forward Compton 
scattering off a hadron target is represented through two blobs connected by the two-parton intermediate states 
involving either quarks or gluons. Strictly speaking, the convolutions in Eq. (j4|) do not correspond to the conventional 
concept of the QCD factorization because each of the amplitudes there (each blob in Fig. [Ij can contain the both 
non-perturbative and perturbative contributions, so we called them the primary convolutions and considered their 
integrability in the Born approximation. Accounting for the radiative corrections converted the Born amplitudes into 
Eqs. (|33l34p . These convolutions complemented by the restrictions in Eqs. ([551155)) rigorously correspond to the QCD 
factorization concept: the upper blob in Fig. [1] becomes free of non-perturbative contributions while the lower blob 
is non-perturbative. However, the factorization convolutions in Eqs. (|33I34[) differ from the convolutions in the kx- 
and collinear factorization and involve the totally unintegrated parton distributions, so we call such a factorization 
the basic factorization. The convolutions in Eqs. p3l34p are UV-stable when the restrictions (|41l43p are fulfilled. 
Applying the Optical theorem to Eqs. p4l33p . we arrived at the expressions (|37l38p for the singlet and non-singlet 
structure functions in the basic (totally unintegrated) form. 

We showed that reducing the basic factorization to the kx- factorization can be done only when the approximation 
in Eq. (|44p was accepted. Using this result, we obtained Eqs. (|50l55p for the structure functions in the framework 
of the kx -factorization. Investigating the impact of the integrability restrictions (|41I43I) on (|50l55p . we obtained the 
theoretical restrictions on the parton distributions at the the kx -factorization: Eqs. (J59l60p show that the non-singlet 
parton distributions should be regular in the longitudinal momentum /3 while the distributions for the singlet Fi can 
include the singular factors /3~'', but with h < 1. Besides, the parton distributions should be regular in the transverse 
nromenta. These restrictions should be respected in fits for the parton distributions composed in the framework of 
the kx -factorization. When the fc_L-dependence of those parton distributions exhibits a sharp maximum at k'j_ = /i^, 
as shown in Figs. 5,6, the kx -factorization can be reduced to the collinear factorization, with fj, playing the role of 
the factorization scale. This our treatment of /i can be checked by fitting experimental data in the framework of the 
kx- factorization. 

The integrability of the forward Compton scattering amplitudes at the basic factorization leads to the restrictions 
(I64I65P on the singular factors x^"- in the fits for the initial parton distributions: a < 1 in the fits for the singlet Fi 
and a ^ in the fits for other structure functions providing the perturbative corrections are attributed to the upper 
blob. 

There is no essential difference between the restrictions in Eqs. (|64I65|) and (|59|60p obtained for the collinear and 
kx- factorizations respectively. On the other hand, the restrictions (|60I65I) on the unpolarized singlet differ from the 
restrictions (|59l64p for all other parton distributions. However, this difference results altogether from the peculiar 
(~ ^/x) perturbative x-evolution of Fi singlet compared to the other structure functions, i.e. from the difference 
between singlet and non-singlet coefficient functions. So, as the same parton distributions contribute to the singlet 
Fi and the other spin- independent structure functions, all the parton distributions should satisfy Eqs. (|59l64p . The 
same is true for the polarized parton distributions. The necessity to keep the singular factors in the fits should be 
regarded as a clear indication that the perturbative components of the structure functions used in the analysis lack 
the resummation of the contributions ^ In (l/x). 

To conclude, let us notice that the results for the non-perturbative contributions to the Compton amplitude and 
DIS structure functions obtained in the present paper follow from general properties of the primary convolutions in 
Eq. Q and because of that they can easily be applied to other inclusive hadronic processes at high energies. 
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Appendix A: Simplification of 5^,^ 

The tensor kAjj',y k in Eq. (|10p depends on the Lorentz subscripts /i and f through the term 

B^^ = k"f^{q + k)"ff,k. (Al) 

-B^^ contains contributions symmetrical and antisymmetrical with respect /i and u: B^^ = B^f,^" + B^J^"^ . Let us 
consider first the antisymmetrical part, Bj^J^"^ . It contributes to the spin-dependent part Aj^^^^ of the amplitude 
A]^^^ . Bfnf^"^' can be simplified as follows: 

Qiasym) ^ _,(! + a)e^.Xpqx [75^7pfc] • (A2) 

The symmetrical part of i?^^ can be simplified similarly. For example, the contribution in B^f^"^ proportional to g^^ 
is 

g^^(l + a)qp[k"fpk]. (A3) 

Now let us simplify kjpk. 

Using the Sudakov parametrization k = aq + /3p + k± in the c.m.f. where p = y/w/2{l, 0, 0, 1) and q = 
\/w/2{\, 0, 0, —1), we rewrite k^pk as follows: 

ki_ipk^ = (fc±7pfc±)^ + {ki_ipki_f, (A4) 

{k±jpk±)^ == (ag + I3p)-ip{aq + /3p), 
{k±jpk±)'^ = k±jpk± 

As the subscript p runs in the longitudinal space (p = 0, 3), 

(kjpkf = A:i7,. (A5) 

This is the only important term in LLA. When we go beyond the LLA, the longitudinal term cam also be essential. 
In order to deal with it, we represent jp in terms of q,p: 

p+q p-q 

7o = —^, 73 = -^=^ (A6j 



Substituting it into Eq. (IA4[) and neglecting terms ^ q'^,p'^, we obtain 



a^qpq + P'^pqp _ ^^^^^^^ ^ ^2^„, , / „2 
/w 



ik^j^k^)^ ^ '^^irL_^j:2^ ^ ^[(^^ + p^)^^ + (_«^ + ;3^)^3], (A7) 



(1 t \L a^qpq-P'^pqp ,,2 o2n , , 2 , r.2^ 1 



The contributions ~ 0^ in Eq. (|A7p can be dropped when Leading Logarithmic Approximation is used to calculate 
^(pe''*)_ It follows form Eqs. (|A41 IA7p that k^pk ~ a^ when a — ?► 00. On the contrary, when a is small, k^pk = k\^p 
. It is convenient to use the interpolation expression for k^pk: 

hpk-lp[w{a^+l3^) + kl]=jpB. (A8) 

It leads to the following approximative expressions for the symmetrical and antisymmetrical parts of Bp^^: 



Blj"'^ = g,.qB, (A9) 

4-^") = ^ep,xpqxlpB 
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Appendix B: Projection operators for the Compton amplitude 



The standard way of studying Afj_i, is to introduce the projection operators and expanding yl^^ into a series of 



invariant amphtudes. The spin-independent part Aj^"^" of A^^ is parameterized as follows: 



where 



A{unpol) ^pil)A + P(2)4„ 



(Bl) 



p(i) 



-.9p 



, M^ ] p{2) 



p _ TO 



P -a^ 

9 



1 

pq 



(B2) 



while the expansion of the spin-dependent part A^ 



{spin) 



IS 



j^(spin) ^ p(.S)A ^ p(f}A. 



(B3) 



with 



-''Vijy — ^^fifXpQX'Jp 1 ^au — ^^fiiyXpQX'^p 



pq 



pq 



^p Pp 



qS 
pq 



(B4) 



Any of the Compton invariant amplitudes ^1,2,3,4 includes both the flavor singlet and non-singlet contributions. Their 
imaginary parts correspond to the singlet and non-singlet components of the DIS structure functions. 
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